Weak- to strong pinning crossover 
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Material defects in hard type II superconductors pin the flux lines and thus establish the 
dissipation-free current transport in the presence of a finite magnetic field. Depending on the density 
and pinning force of the defects and the vortex density, pinning is either weak-collective or strong. 
We analyze the weak- to strong pinning crossover of vortex matter in disordered superconductors 
and discuss the peak effect appearing naturally in this context. 



Pinning of vortices by material defects is crucial in es- 
tablishing the defining property of a superconductor, its 
ability to transport electrical current without dissipation. 
Collective pinning theory [lj, describing the concerted 
action of many weak pins on the vortex system, is play- 
ing a central role in our understanding of this complex 
statistical mechanics problem 0. On the other hand, 
first attempts describing flux pinning go back to Labusch 
0, who described the interaction between vortices and 
strong pinning centers which introduce large (plastic) de- 
formations in the vortex system. In this letter, we de- 
scribe how these two theories relate to one another; given 
the density n p and force / p of pinning centers, as well as 
the vortex density n v = 1/oq, we identify the regimes 
where individual vortex lines and the bulk vortex lattice 
are pinned by the collective action of many weak pins 
or by the independent action of strong pins, see Fig. ^ 
We naturally recover the peak effect |j] described in the 
work of Larkin and Ovchinnikov Q and establish its for- 
mal relation to the Landau theory of phase transitions. 




In a type II superconductor, the field (B) induced vor- 
tices subject to a current flow j experience the Lorentz 
force density F L = j AB/c and the resulting vortex mo- 
tion leads to dissipation. The superconducting response 
is resurrected through material inhomogencitics pinning 
the vortices at energetically favorable locations. The pin- 
ning force density F p ; n defines a critical current density 
j c = cFpi n / B below which the current can flow free of dis- 
sipation. Usually, this critical current density is consid- 
erably reduced with respect to the depairing current den- 
sity j ~ cH c /Att\ ~ ce /$o£; here, H c = $ /2\/27rA£ 
is the critical magnetic field, A and £ denote the pene- 
tration depth and the coherence length £, respectively, 
$0 = hc/2e is the flux unit, and £q = (^q/AttX) 2 is the 
(line) energy scale. Below, we focus on the most generic 
situation of isotropic superconductors and ignore effects 
due to thermal fluctuations. 

When pinning is strong (HEI-Q defects act individually 
and the pinning force density F p is linear in the density 
n p and average pinning force (/ P i n ) of defects. The clas- 
sic arguments characterizing strong pinning go back to 
Labusch 0] , see also 0, ®] : A strong pinning defect in- 
duces plastic deformations in the vortex lattice 

BUS) 

and the energy landscape e P i n (R) becomes multi- valued 



FIG. 1: Pinning diagram delineating the various pinning 
regimes involving collective- versus individual pinning and 
ID-line- versus 3D-bulk pinning (/Lab denotes the Labusch 
force): 3D wep - bulk weak collective pinning, ID cp - col- 
lective line pinning, ID sp - strong line pinning, 3D sp - bulk 
strong pinning. The insets illustrate the weak- and strong 
(plastic) distortions of the vortex system characterizing the 
pinning regimes. Lines refer to crossovers. 



in the vortex position R, see Fig. El The averaging over 
defect locations then has to account for the preparation 
of the system. We concentrate on the critical current 
density and thus search for the force against drag; the 
vortex position then is parametrized through the two- 
component drag parameter Rd fixing the position of the 
unperturbed lattice with respect to the defect. Dragging 
the system along the x-direction, we express the drag 
force — d x e p i n (x, y) integrated along x through the jump 
Ae p i n (y) > in the pinning energy and average over 'im- 
pact parameters' y, 



(/pin) — 



dx 



Jy dy d x e pin (x,y) = _ p A^y) 
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where a denotes the distance between periodic jumps |9(. 
For moderately strong pins with deformations not ex- 
ceeding the lattice constant we have a sa ao and assum- 
ing a maximal transverse trapping distance t± along the 
y-axis we obtain the mean pinning force 
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FIG. 2: Energy landscape e p i n and pinning force / p i n versus 
displacement Rd of the vortex lattice relative to the defect; for 
weak pinning these are single- valued functions in (dashed 
lines) , while strong pinning produces plastic deformations and 
renders e plrL , / p ; n multi-valued (solid lines; dotted lines indi- 
cate unstable branches). Bottom right: Trapping geometry 
(top view) for a circularly symmetric situation. 

with the jump Ae p i n (0) ~ ty/ p expressed via the typical 
impurity force / p and the bistability range tu of e p i n (x, 0); 
the product fj_t|i defines the trapping area <St rap associ- 
ated with the strong pin Q. The low impurity concen- 
tration n p implies non-interfering defects and we obtain 
a critical current density j c — — cn p (f p i n ) / B linear in n p , 

j c fa (c/B)n cB f p fa j [n p £S tIap ] f p /e , (2) 

with the effective impurity density n e g — n p (St Tap / ci^) . 

In order to derive a quantitative criterion for the ap- 
pearance of strong pinning, we consider a single defect at 
the origin with a pinning potential e p (r). Such a defect 
acts on the vortex system to produce a pinning energy 
density E p (r,u) = J2„ e p (r)J 2 (R — R y — u(R„, z)), with 
vortices positioned at R^ + u(R„, z), H v the equilibrium 
positions and u the displacement field. The latter follows 
from the solution of the implicit equation (r v = (R y , z)) 

u a (v v ) = J G?r'G a0 (Y v -v r )[-d^E p ](r\u') 
= J2[ dz>G aP (r„ - r' v )fpAK + u(ii), z>) 

v' 

= G a p(R v - R d , 0)/ p>/3 (R d + u(R d , 0), 0), (3) 

with G a p(r) the elastic Green's function and f p = 
— V„e p (u) the pinning force of the defect. In the last 
equation we have assumed a moderately strong pinning 
potential (pinning one vortex at most) of range much 
smaller than the lattice constant ao and have chosen R d 
as the distance to the vortex closest to the defect [l^. 
Evaluating @ for r v — (R d ,0), we arrive at the self- 
consistency equation (note that G a p(r = 0) is diagonal) 

u(R,0)«C- 1 f p (R + u(R,0),0), (4) 



with the effective elastic constant C 1-1 = J d 3 fc/(27r) 3 
G xx (\i). For weak pinning the displacement u is 
small and the solution u(R, 0) « f p (R)/(7 is unique. 
Strong pinning, however, produces multi-valued func- 
tions u(R,0) and e pin (R), cf. Fig. [3 The solution of Q 
turns multi-valued as the displacement collapses when 
diju — > co. Assuming a defect symmetric in the plane, 
e p (R, z) = e p (R,z), and dragging the lattice through 
the defect center along the ir-axis, we find u' — f p (x + u) 
[C — f p (x + u)]^ 1 (note that x > implies u < 0) and 
arrive at the (Labusch) criterion in the form 

d x f P - -8 2 x e p = C; (5) 

hence, in order to produce strong pinning the (negative) 
curvature of the pinning energy e p has to overcompensate 
the lattice elasticity (the Labusch criterion involves the 
maximal negative curvature above the inflection point). 
Note that the Labusch criterion tests an individual pin- 
ning center and classifies it as a weak or strong one. 

When pinning is weak, the elastic forces dominate over 
the pinning forces and the defects compete; we then are 
faced with the problem of the statistical summation of in- 
dividual pinning forces. For weak pins the average (/ P i n ) 
vanishes and pinning is due to fluctuations in the pinning 
force density: the forces of the competing pins (with pin- 
ning force f p , density n p , and extension r p ~ £) add up 
randomly and produce the pinning energy 

K^V)) 1/2 ^[f 2 p n p ^/a ) 2 V] 1/2 i; (6) 

only vortex cores are pinned by the disorder, hence the 
factor (£/ao) 2 . Within weak collective pinning theory 
the sublinear growth of (E^CV)) 1 ^ 2 with volume turns 
linear when the displacement u increases beyond the scale 
£ of the pinning potential, thus defining the collective 
pinning volume V c - Each volume of size V c is pinned 
independently with a pinning energy U c — (Sp^Vc)) 1 ' 2 
and we obtain a proper pinning force density 

F p - m ~ U c /V c r p ~ (f 2 p n p (i/a a ) 2 /V L ) 1/2 - (7) 

balancing this pinning force density against the Lorentz 
force density jB/c we find a finite critical current density 
j c ~ aFpin/ B. The remaining task is the determination 
of the collective pinning volume V c ; its calculation is com- 
plicated by the dispersion and anisotropy of the vortex 
lattice, see below and Ref. [2] for a detailed discussion. 

It is instructive to compare the weak- and strong pin- 
ning schemes and their dependence on dimensionality, 
particularly in the limit of a small defect density n p (in 
the following, we assume pinning sites characterized by 
their force / p and extension £). An isolated vortex line 
(ID) is always subject to strong pinning forces as the ef- 
fective elastic coefficient C vanishes due to the diverging 
integral. At the same time, the deformation of the line 
due to the pins is large and we cannot ignore their mu- 
tual competition. Comparing the elastic energy e^ 2 / L c 
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and the pinning energy U c = (fp-n-pLc^ 2 ) 1 / 2 ^, we find the 

collective pinning length L c ~ (e"o/ f P n pY^ 3 an d a critical 
current density 



ic-joKr/p/^o) 



,2/3 



(8) 



This result is valid as long as many pins compete within 
the volume £ 2 £ c ; the condition n p ^ 2 L c > 1 defines the 
lower limit n 1D ~ / P /£o£ 3 < n p where the critical current 
density assumes the value j c ~ Jo(/p/eo) 2 - 

For small densities n p < n 1D the pins acts individually 
and we determine j c from the force balance (&o/c)j c lu ~ 
Ae P i n ~ fpU, with u ~ ty the displacement directed 
along the force. The displacement u and the length I 
between two subsequent pins fixing the vortex derives 
from an analysis of the pinned vortex geometry, see Fig. 
^ inset: integrating the force equation £q8 2 u = f(z) 
(with f(z) the force per unit length acting on the line) 
over one pinning center, we find the distortion angle 
9 = d z u ~ u/l ~ fp/eo 10] . A vortex segment of length 
I deformed by the angle 9 in the direction of the driving 
force encounters 9l 2 £n p defects (with the trapping length 
t± ~ £)• At the distance I, this number is unity, hence 
I ~ -\/ e o//p n p^ ana - we obtain the critical current density 



Jc ~ joKC 3 /p/eo) 



»V2 



(9) 



At the crossover density n 1D ~ / P /eo£ 3 the critical 
current density matches up with the weak pinning re- 
sult; also, the displacement u ~ I fp/eo is of order £ at 
the crossover density n 1D and hence matches the dis- 
placement field relevant in the collective pinning sce- 
nario. Note that collective pinning (JSJ dominates over 
the strong pinning J^J at large densities n p > n 1D . 

For the vortex lattice (3D bulk pinning; we assume 
ao < A) the Labusch criterion J^J offers a distinction be- 
tween weak and strong pinning centers; using the Green's 
function for the vortex lattice (see, e.g., 0) we find 
C ~ eo/ag. According to (J5J a pinning center changes 
from weak to strong at f p ~ /Lab = £oC/ a o- We first 
review the weak pinning situation with f p < /Lab (where 
necessary, we encode quantities in this regime with a su- 
perscript '<'). The determination of the anisotropic col- 
lective pinning volume V c = R 2 C L\ has to account for 
the dispersion in the tilt modulus at intermediate scales 
a < R c < A, see Ref. Q, and produces the results 



. . e 
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Jc ~ JO ^ 



«0 

L c 
ao 

Lr 



e -2c[L c /a ] 3 i Rc < Aj (1Q) 



Rc > A; 



(11) 



we have made use of the single vortex pinning parame- 
ter Lc/ao ~ (/^ab/ ' fpao^np) 1 ^ 3 '■ The numericals c and 
v follow from a 2-loop renormalization group analysis 
[ill Xv\ . These results are valid as long as many (com- 
peting) pins are active in the volume V c , n p (£, 2 / Qq)V c > 1. 



For f p < /L a b this condition is violated in the large 
n p limit. However, with increasing pinning density n p , 
the collective pinning radius R c decreases, first below 
A at n x ~ flab/ fp a o£ 2 m (A/a ) delineating the disper- 
sive regime, and then below ao at n< Q ~ flab/ fp a o£, 2 
where the condition n p (£ 2 / 'a^)V c > 1 is still fulfilled. At 
the crossover density the 3D weak collective pinning 
crosses over to the ID weak collective pinning result (jHJ. 

Turning to strong pinning / p > /L a b (encoded with a 
superscript '>') we start at high densities; as the Labusch 
criterion is not effective in ID, the system remains collec- 
tively pinned for n p > n lrt and crosses over to ID strong 
pinning as n p drops below 

Tim- Decreasing tl p further, 
the pinning distance / ~ ao yj (/Lab / fp)/ ^p«o£ 2 increases 
beyond a as n p decreases below n> ~ (/Lab// P )/ao£ 2 
and the system enters the 3D strong pinning regime, see 
Fig. ^ The calculation of the mean pinning force den- 
sity F p i n ~ n p (f p i n ) proceeds along the lines discussed 
above and involves the trapping area Strap ~ with 
t± ~ £ and t|| ~ u ~ fp/C; we obtain the force density 
Fp- m ~ n p(C/ a o)/p/ £ an d a critical current density 



f 2 

Joao£ n p — 



f 2 

jo -o « P aoC 



(12) 



Lab 



The bulk strong pinning result ({12(1 smoothly transforms 
into the ID expression at n> where I ~ ao. On 
the contrary, the strong pinning expression 112|) does not 
match up with the bulk weak collective pinning results 
(fT7Tf> and (fTT|l at / p = /L a b (we concentrate on low im- 
purity densities with n p ao£ 2 < 1, cf. Fig. However, 
we have to keep in mind that our rough derivation of the 
strong pinning result l|12|) breaks down on approaching 
the critical force /Lab- Indeed, since the displacement 
field u(r) turns single valued below /Lab, strong pinning 
vanishes altogether (with a power [/ p — /Lab] 2 , see |T5]l) 
and pinning survives only in the form of weak collec- 
tive pinning due to fluctuations in the impurity density. 
Within the approximative scheme adopted here the sharp 
rise of the critical current density at f p > /Lab is encoded 
in a jump j c | sp /j c | W cp ~ (A 2 /a§)/n p a ^ 2 > 1 for n p < n x 
(~ exp[2c/n p ao£ 2 ] for n p > n\). 

The crossover from strong to weak pinning at the 
Labusch condition JSJ can be analyzed within a Lan- 
dau type expansion: We define the free energy functional 
e p j n (u, Rd) = Cu 2 /2 + e p (Rd + u) from which the self- 
consistency equation (jJJ follows by variation. Note that 
the derivative — d x e pm = f p>x (R.d + u) provides the force 
along x acting on a vortex separated from the defect by 
Rd and deformed by u, cf. Fig. 1; it is this force which 
has to be averaged over in the definition of (/ p i n ). 

We first concentrate on the trajectory Rd = (x, 0) with 
u = (u, 0). The curvature e''(u) relevant in © assumes 
a maximal negative value; we denote the corresponding 
location and value by u K and —k, respectively. Next, we 
expand the curvature around u K , e p (u) ~ — k + a(u — 
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u K ) 2 /2; integrating in u and combining with the elastic 
term Cu 2 /2 we arrive at the expansion 

e p i n [w, x] ~ Cu 2 /2 + v (x + u — u K ) (13) 
- k(x + u - u K ) 2 /2 + a(x + u - u K ) 4 /24. 

This pinning energy maps to the free energy e mag [<^, h] = 
T(j) 2 /2 + a<^> 4 /24 — hcf> of a one-component magnet in 
a magnetic field 01 if we define the 'order parameter' 
<t> — x + u — u K , the 'temperature difference' r = C — k, 
and the 'magnetic field' h = C(x — u K — v/C). The 
'high-temperature' phase r > describing the param- 
agnet corresponds to weak pinning, while the two fer- 
romagnetic phases at 'low temperatures' r < stand 
for the pinned (</> < 0) and unpinned (</> > 0) states; 
the transition between these states is discontinuous and 
the associated coexistence regime extends over the 'field' 
domain \h\ < h* = (2/3(7) ^/2/^\t\ 3/2 . At h* the trap- 
ping/detrapping of the vortex from the defect produces 
jumps A</> = 3^2/a\r\ 1 / 2 in the 'order parameter', lead- 
ing to jumps Ae = Ae p i n /2 = (9 /2a) t 2 in the energy. 

For finite 'impact parameters' y we have to determine 
the trapping distance t±. Assuming rotational symme- 
try, the bistable regime is bounded by a circle of radius 
R = x* = u K + v jC and hence t± w 2x* (note that at 
r = we have h* — but the critical drag parameter x* 
does not vanish). The (uncompensated) trapping area 
determining the average pinning force (/pin) is shown in 
Fig. [2 combining the above results for the jump in pin- 
ning energy and the transverse trapping distance we find 
the averaged pinning force (c.f. (JIJ) 

(/pin) « 18K + v/0)[C - n] 2 /aa 2 . (14) 

Defining the individual force of (equal) pinning centers 
via / p = max u [d u f p ](u)£ = k£ (then / Lab = Cg) we can 
translate (I14|) into an expression for the critical current 
density j c extending the strong pinning result 112fl to the 
vicinity of the Labusch point, 

3c ~ Jo (e/a 2 Q )n p a Q e [/ P // La b - l] 2 - (15) 

Comparing with the weak pinning result (II 1|) . we note 
a sharp rise in the critical current density j c once the 
strong pinning force overcomes the weak pinning result 
0. With the small parameter S — (a /X)y / n p aQ^ 2 < 1, 
this crossover appears above but still close to the Labusch 
point as /Lab oc C decreases below / p /(l + S). 

Another remarkable result is the interpretation of the 
collective pinning scenario in terms of the strong pinning 
picture; indeed, summing over competing pins within the 
collective pinning volume V c produces the correspond- 
ing critical Labusch force. Quantitatively, we compare 
the force gradient /' ~ [n P (£, 2 /afyV] 1 / 2 f p /t; accumulated 
within the (anisotropic) volume V = LR 2 = (A/ao)i? 3 
with the elastic parameter C(R) — SoAR/clq for smooth 
distortions on the scale R > A (non-dispersive regime) 



and apply the Labusch criterion (0. We then find the 
scale R c = Xfl^/ f p n p ao^ 2 , where the accumulated pin- 
ning force overcompensates the elastic force; this length 
agrees with the 3D collective pinning length in the non- 
dispersive regime Q ■ The resulting bistable solutions are 
the signature of the alternative pinning valleys which the 
collective pinning volume can select beyond the scale R c . 

The above discussion sheds light on the general concept 
of pinning. Pinning is absent in the rigid limit. A finite 
but large elasticity (with /Lab > /p) allows only for weak 
deformations and individual pins cannot hold the lattice 
as the averaging over individual pinning forces produces 
a null result. Hence, pinning is due only to fluctuations 
in the pinning forces and thus collective. Reducing the 
elasticity, strong pinning defects appear when /Lab drops 
below / p ; they pin the lattice individually and strong 
pinning, linear in the defect density n p , outperforms col- 
lective pinning. The important role played by the cur- 
vature e p < in the pinning potential is an interesting 
topic for numerical studies. The crossover between weak 
collective and strong pinning can be realized in experi- 
ments: increasing the magnetic field towards its critical 
value H C2 leads to a marked softening of the elastic mod- 
uli. The reduction in the elastic moduli entails a decrease 
of the Labusch force /L a b an d triggers the crossover from 
weak- to strong pinning, producing the well known peak 
effect in the critical current density Q, 0] . 
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